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1. Introduction 

The dynamics of relativistic thin shells is a recurrent topic in the literature about the 
classical theory of gravitating systems and the still ongoing attempts to obtain a coherent 
description of their quantum behaviour. Certainly, a good reason to make this system 
a preferred one for a lot of models is the clear, synthetic description of its dynamics 
in terms of Israel's junction conditions [I] [2] (the null case, considered in detail in the 
seminal paper of C. Barrabes and W. Israel [Sj, is also interesting for null-like surfaces 
[3], i.e. light-like matter shells [S], for which an Hamiltonian treatment is given in [H] 
generalizing the approach described in Using this formalism!, which has an intuitive 
geometric meaning, many relevant aspects of gravitation have been brought into light. 

Gravitating shells have indeed been considered as natural models for different 
astrophysical problems: the description of variable cosmic objects [0] and of specific 
aspects (like ejection [TU] or crossing of layers [TTj, critical phenomena [TJ], perturbations 
[T3] and back- reaction ^3]) in gravitational collapse ^3103 E] are only a few examples. 

Moreover, at larger scales, specific configurations of shells have also been considered 
to construct cosmological models lB i (even with hierarchical (fractal) structure [19J), 
to analyze phase transitions in the early universe j2H| or to describe cosmological voids 
[23] ; semiclassical models have tackled the problem of avoiding the initial singularity of 
the Big-Bang scenario by quantum tunnelling [221 E3 E3j • 

As a matter of fact quantum semiclassical models have conveniently been employed 
as useful simple examples to better understand possible properties and modifications 
to the spacetime structure at scales at which quantum effects should give significant 
contributions to gravitational physics. Apart from the quantization of the gravitating 
shell itself (as in [213 EE]), considered also in the context of gravitational collapse |27j . 
models have been proposed to study quantum properties of black holes [23 EH EH] and 
their formation process [HI] as well as to analyze wormhole spacetimes [H21 EH1 E31 EE] 
and the quantum stabilization of their instability [HEIEU], targeting the fuzzy properties 
of spacetime foam [HE] and Planck scale physics [HH] • 

Other problems of fundamental nature in quantum gravity have received attention 
through the study of shell dynamics: as an exemplificative list, we mention here Hawking 
radiation j301EniE] the horizon problem in wormhole spacetimes [32], the time problem 
in canonical relativity [43J, the problem of localization of gravitational energy [33], the 
thermodynamics of self-gravitating systems ^El HD) and the possibility of connecting 
compact with non-compact dimensions [TF] . 

Many of the above discussions have been performed under the simplifying assumption 
of spherical symmetry: this is especially useful in the quantum treatment, because 
the minisuperspace approximation greatly reduces the complexity of the mathematical 
treatment. 

But, at least at the classical level, studies have also been performed for cylindrical 

| But see also [H] and references therein for a complementary approach which also tackles the issue of 
stability. 
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models (see e.g. [3E1 HE ESQ ) • 

With the development of the models shortly cited above, particularly those 
involved with the quantization of the system, many subtleties emerged as byproducts 
of corresponding difficulties already encountered in tentative approaches to Quantum 
Gravity and mainly related to the reparametrization invariance of the theory. Since the 
junction conditions, essentially, are a first integral of the equations of motion of the shell, 
many authors revolved their attention to the derivation of these equations starting from 
an action principle. A consistent Lagrangian/Hamiltonian formalism has been developed 
jSIl E3 LH3 E3] (also reduced by spherical symmetry [53]), and the relevant degrees of 
freedom of the system [HE] discussed together with a variational principle, which is also 
the subject of j^ZlEHl (interesting considerations can also be found in [HH|). 

Recently even more interest in the thin shell formalism is coming thanks to the 
development of brane world scenarios, where our universe is seen as a four dimensional 
brane embedded in a five dimensional space [60, 61j . This configuration can be given 
a wormhole interpretation [E2] and has also been analyzed from the point of view of 
energy conditions [EB] {not) satisfied in the higher dimensional background. 

In these and other studies, different cases of junctions between spacetimes have 
been considered: for example between anti de Sitter and anti de Sitter [E2], Friedmann- 
Robert son- Walker and Friedmann- Robertson- Walker [33], Minkowski and Minkowski 
[HZ1 EH], Schwarzschild and Schwarzschild [HE1 HU EH IE] , ReiBner- Nordstrom and Reifiner- 
Nordstrom [TSJ ESJ ? de Sitter and Reifiner-Nordstrom [E3], de Sitter and Schwarzschild 
[E31 ESI EI], de Sitter and Schwarzschild- de Sitter de Sitter and Vaidya [24J, 

Friedmann-like and ReiBner- Nordstrom [EE], Minkowski and Friedmann [2J, Minkowski 
and ReiBner-Nordstrdm [131311122!, Minkowski and Schwarzschild |2Hl EH EHl IHl 113 EH! , 
Minkowski and Vaidya [T2], Schwarzschild and Reifiner-Nordstrom [67J, Schwarzschild 
and Schwarzschild-anti de Sitter [3Z], Schwarzschild and Vaidya [3E], Tolman and 
Friedman |19j . Lemaitre-Tolman-Bondi and Lemaitre-Tolman-Bondi [TH] . 

In this paper we are also going to use a general relativistic shell to analyze, even 
if only at the semiclassical level, the problem of quantization of a gravitational system. 
We will restrict ourselves, as it has been done in many of the papers cited above, to 
the spherically symmetric case and we will study the semiclassical quantum dynamics 
in the case in which the shell separates an interior spacetime of the de Sitter geometry, 
from an exterior of the extremal Reifiner-Nordstrom type. An observer crossing the 
shell will naively see some non-vanishing vacuum energy density to be converted into 
physical properties like charge and mass. From the classical dynamics there are no 
restrictions on the values of the physical parameters characterizing the geometry of 
spacetime. But, starting from a Hamiltonian description of the shell dynamics, we will 
try to analyze its quantum behaviour. Lacking a full theory of quantum gravity, which 
would of course be the natural setting for this kind of problem, we will tackle it only 
at the semiclassical level: under this word, we will understand that the action for the 
shell is given as an integer multiple of the quantum, h. We will see that this condition 
results in a constraint on the parameters for the interior and exterior geometries. This 
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is hardly surprising: indeed a full quantum theory of gravity, would have the task of 
determining the probability amplitude for a given configuration of the three-geometries 
taught as points in superspace; in our quantum minisuperspace approach, the only 
free parameters remain the constants (de Sitter cosmological horizon, charge and mass) 
fixing the interior and exterior metrics, and is thus as a relation among them that the 
semiclassical quantization conditions realizes itself. 

With the above ideas in mind the paper is organized as follows. 
In section 121 we will set up our model by giving all relevant definitions; we will also recall 
some well known results adapted to our special case, to fix notations and conventions, 
and will present all relevant dynamical quantities for the computations that follows. 
The Bohr-Sommerfeld quantization condition is also recalled. Then, in section |3] the 
classical dynamics of the system is sketched and the associated spacetime structure 
discussed, with particular emphasis on the bounded trajectories. This prepares the 
ground for section HI where the classical action is numerically evaluated for bounded 
trajectories. This result is then used in section to show how the Bohr-Sommerfeld 
quantization condition characterizes the properties of the semiclassical quantum system. 
After a preliminary rough estimate (subsection EH}, we present the semi-classical results 
for the quantum levels of the shell and the corresponding internal/external geometries 
and approximate the results with a properly chosen analytic (polynomial) expression. 
Discussion about the results and possible refinements of the model follow in section |H1 
Five short appendices are devoted to a more detailed analysis of some technical points. 
The turning points of the classical motion are discussed in |Appendix A| The issue 
about the stability of the classical solution against single particle decay is studied 
in |Appendix B| |Appendix C shows that the bounded trajectories are not affected by 
change of direction of the normal to the shell trajectory. The characterization of the 
singularity that appears in the integral for the computation of the classical action as an 
integrable one is done in |Appendix D and the determination of the leading terms in the 



integrand of the same computation is the topic of Appendix E 



2. Preliminaries 



In this section we define the system, motivate the settings under which we study its 
classical and semiclassical dynamics and recall some useful results and definitions. 

Let us thus start with the geometrodynamical framework, by considering two 
spacetime domains joined along a spherically symmetric timelike shell. We assume, 
for the region we shall call the interior, a geometry of the de Sitter type |6SJ IHHJ EH] (we 
denote with H the cosmological horizon), so that the metric in static coordinates is: 

9in = dia g (/in(r),/ i ; 1 (r),r 2 ,r 2 sin^) (1) 
/in(r) = 1 - j^. 

For the exterior region we choose a spacetime of the Reifiner-Nordstrom type |7T] 1721 170]. 
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with metric given by 

gZ = diag (/ out (r), /"», r 2 , r 2 sin 5) (2) 

/out(r) = 1 — + — , 

M being the Schwarzschild mass and Q the electric charge. Furthermore we join the 
two regions along the timelike trajectory of a spherical dust shell of constant total 
mass-energy m. 

As is well known jH |2j the dynamics of the compound gravitational system is 
encoded in Israel's junction conditions: they match the jump in the extrinsic curvature 
due to the different spacetime geometries on the two sides of the shell surface and the 
(singular) stress-energy tensor of the shell itself. Under the simplifying assumption of 
spherical symmetry considered here, it is possible to reduce them to the single scalar 
equation E3j 

M = ^ (3) 

where as customary we use square brackets as a shorthand for the jump of the enclosed 
quantity in the passage from the "in" to the "out" domain across the shell§, i.e. 



[X] :— X in — X, 



out; 



and 



R 2 

( a P)in := ff inAn = a in \f R 2 + 1 - — 



2 



Mout : = ^outAmt = 0- out \jR 2 + 1 - ^ + ^ . 

In the above expressions R = R(t) is the shell radius expressed as a function of the 
proper time r of an observer co-moving with the shell and we denote with an over-dot 
the (total) derivative with respect to r. 

From the general theory of shell dynamics, we know that the signs of the radicals 
do matter, being related both to the side of the maximally extended diagram for the 
spacetime manifold, which is crossed by the trajectory of the shell [H3], and to the 
direction of the outward (i.e toward increasing radius) pointing normal to the shell 
surface [HI2|- This is the reason why they are denoted explicitly by o" in / ou t- Their values 
can be analytically determined thanks to the results [Z31 HI] 

= a ia (R) = -Sign jm ^ - 2MR + Q 2 - m 2 ) j (4) 



1 



o-ont = o- out (R) = -Sign jm (J^ - 2MR + Q 2 + m 2 J j , (5) 

which can be obtained by properly squaring the junction condition Q. 

§ To avoid any possible confusion, in what follows we are going to use square brackets only with this 
meaning, according to the following definition. 
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Following the notation of reference [73] we know that the junction condition (J3J) can 
be derived as the Superhamiltonian constraint for the system, where the corresponding 
Superhamiltonian is then nothing but 

H = PR- C = R - m, (6) 

L being the Lagrangian density 



m 



R 



a/3 - -R In 

H 2 



ap + R 



a(3-R 

and P being the conjugate momentum to the canonical variable R: 

a(3 + R 



dR 



In 



(7) 



(8) 



a(3-R 

The dynamics of the system can be studied with the help of an effective equation of 
motion, which is useful in removing the square roots in (J3J) and can be put in the form 
of a classical one-dimensional dynamical problem, the motion of an effective particle of 
unitary mass and vanishing total energy |M] 175] 

R 2 + V(R) = (9) 

in a potential given by 

R 8 - 4H 2 MR 5 + 2H 2 {m 2 + Q 2 )R 4 



V(R) 



4m 2 H 2 R 2 
4H 4 (M 2 - m 2 )R 2 + 4H 4 M(m 2 



Q 2 )R + H 4 (m 2 - Q 



2\2 



4m 2 H 2 R 2 ^ 
To evaluate the classical action along a classically allowed trajectory, we need an 
expression for the effective momentum P evaluated along the same trajectory. Thus 
we have to substitute for the R dependence in (JSJ) and using in this procedure relation 
© we get 



P(R) 



-i?tanh 



with s(r) = Sign 




2mH 2 R^-V(R) 



-R 4 + 2H 2 R 2 - 2H 2 MR + H 2 (Q 2 - m 




(11) 



R?_\ 
H 2 J 



1 



2M Q 2 
~R + R 2 



We can now use the above result to compute the action along a classically allowed 
trajectory, having turning points at Ri and R 2 , since 

P(R)dR, (12) 



classical 



Hi 



and then implement a semiclassical quantization scheme a la Bohr-Sommerfeld, by 
considering allowed quantum states to have the action as an integer multiple of the 
elementary quantum [TH] lp = h = 1: 

^classical = n i 71 = 1,2, ... . (13) 
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To successfully complete this task we need in first place an analysis of the allowed 
bounded classical trajectories, which we will perform in the next section. 

Before embarking this program, let us shortly comment about the semiclassical 
quantization procedure outlined above. There are indeed many different approaches for 
the quantization of gravitational systems, and it is not often clear which should be the 
preferred one. Moreover, deep ideas have already been discussed to a great extent in 
the literature cited above. It is not the goal of this paper to address this fundamental 
problem, but we think it is important to give a short account about the reliability of the 
results that will be derived in what follows. In particular a formalism using expression 
(JTTJ) , but evaluated along a classically forbidden trajectory, has already been successfully 
used in [77] and [73] to reproduce some well known results about vacuum decay and 
the influence of gravity on it, already studied in the seminal papers by Coleman and 
de Luccia [IB] and by Parke [ZH] . Moreover, as already noted by Sommerfeld in the days 
of the early development of Quantum Mechanics [70], the quantization condition (fT3|) 
is 11 particularly valuable, for it could be applied both to relativistic and non-relativistic 
systems" . Thus, we think that the above considerations justify our tentative approach, 
in which the semiclassical quantization condition is applied, through an already tested 
procedure, to a classically well-know gravitational system. 

3. Classical Dynamics 

The results presented above are valid for arbitrary values of the four parameters entering 
the problem, namely the mass M and the charge Q of the external Reifiner-Nordstrom 
spacetime, the de Sitter radius H of the interior geometry and the total mass-energy of 
the dust shell, m. We now specialize them to a more particular setting (which has the 
advantage of removing the second line in expression for the effective potential): 

(i) we take the external Reifiner-Nordstrom spacetime to be extremal, i.e. with 
\Q\=M- 

(ii) we assume that the total mass energy of the shell is m = \Q\. 

The Penrose diagrams for the full de Sitter and extremal Reifiner-Nordstrom [HH] 
spacetimes are shown in figures ^ and [2] respectively. 

Before studying the possible shell trajectories in the two geometries to identify 
the bounded ones, which we are interested in, we take full advantage of the parameter 
reduction implicit in the assumptions above, by passing to adimensional variables: this 
will be more convenient also for the subsequent numerical treatment. We thus choose 
to parametrize all the variables and constants in terms of the de Sitter cosmological 
horizon H by setting 

R r n Q M m 
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i° 



Figure 1. Maximally extended Penrose diagram of the de Sitter spacetime. 




Figure 2. Maximally extended Penrose diagram of the extremal Reifiner- Nordstrom 
spacetime. 
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Then the quantities which are functions of R, become functions of x, all retaining their 
numerical values but P and S, which are rescaled by H and H 2 respectively: 

Hri6)= p&m. , m ^jEM (15) 

We thus have for the quantities evaluated along a classical trajectory, which are relevant 
in the study of the classical dynamics || 

/ in (x) = l-x 2 (16) 
2181 8 2 

U{x) = \- + ^ (17) 
X x z 

( x in (x) = -Sign(x 3 -2|8|) (18) 
°oMt{x) = -Sign(x 4 -2|8|x + 28 2 ) (19) 
4|6|: 
48 2 



? = _s*(*«-4|8j, + 48') 




Sign{(l-x 2 )(a;-|G|) 2 } 



= - Itmh " < i ^rf^r 1 > ; (21) 

of course all can be expressed as functions of the single adimensional parameter 8. 

Following foQ ITS'] we can study the classical dynamics in a compact way by means 
of a comprehensive graphical method fully exploiting the handy relation It consists 
in plotting the potential V(x) together with the metric functions f- m , f out . Then the 
allowed trajectories with the corresponding turning points can be determined looking at 
the segments of the x-axis (corresponding to zero energy), that are above the graph of 
the potential. In this diagram the points where the metric functions vanish, quickly help 
in determining if a classical path crosses the horizons of the external/internal geometry. 
Moreover the rescaled values of the radial coordinate x for which o" in / ou t changes sign are 
given, if they exist, by the values at which the metric function plots are tangent to the 
graph of the potential. This graphical information can be completed by the following 
analytical results. 

Turning points of the potential : 

from (}2T)|) we see that the potential for 6 ^ has a (double) zero at x = 0, so that it 
is regular at the origin, which is thus a trivial turning point of classical trajectories. 
Other turning points may, or may not, be present, depending on the value assumed 
by the parameter 8. Two cases are possible, as shown in figure El Either there can 
be no other turning points, so that only a so called "bounce" classical trajectory 
exists (as is the case in figure El for 8 = 0.7), or there can be two more turning 
points so that in addition to the bounce trajectory there is also a bounded one 
(this is also shown in figure El for 8 = 0.5). As explicitly proved in Appendix A[ 



the critical value for the parameter 6, 8 cr it. = (3/4) 3//2 gives the in between case, 

| We denote with an overbar quantities, let us say g, which are function of the rescaled radial coordinate 
x, although in many cases we have g(x) — g(R). For the sake of precision, note that fi n {R) = fin(x), 
/out = /out(aO) a in(R) = &in(x), cr out (R) = CTout^), and V(R) = V(x). 



WKB states of a dust shell 



10 



V(x) 
0.2 - 


i i i \ i . 












1 ' X 






0.2 


0.4 0.6 0.8 


l.a 


1.2 


-0.2 ■ 












-0.4 - 






= 0.5 | 






-0.6 - 













V(x) 
0.2 - 


=d 1 1 1 1 






0.2 


TW 8*— """"0 1 


X 

1.2 


-0.2 - 










-0.4 ■ 






| = 0.7 1 




-0.6 ■ 











V(x) 






0.2 - 




1 1 =t — ^-i \-~ 




cl2"~ tra — rrefo.8 


m 1.2 


-0.2 - 






-0.4 - 




e = e crit . = (3/4) 3 / 2 




-0.6 ■ 







Figure 3. Graph of the potential V(x) for different values of the parameter 0. 
Depending on the value of 0, the classical trajectory can have two non-vanishing 
turning points, or no non-vanishing turning points, as shown in the first and second 
figure, respectively. The in between case, which occurs for O = cl -it. = (3/4) 3,/2 , is 
depicted in the third diagram. 

when the potential is tangent to the x axis (third plot, again in figure EJ). We 
thus see that only for < < cr it. there are two non-vanishing turning points 
a^min, a^max (actually with x min = x max if 6 = cr it.) an d thus bounded solutions 
are allowed, the classical path being represented by the segment [0, x min ]. We note 
that it is possible to find x m i n and x max in closed form solving the quartic equation 
that gives the non- vanishing solutions of V(x) = 0, and this (not very enlightening) 
expressions are reported in |Appcndix A| 
Horizon positions with respect to the classical path : 

to correctly understand the spacetime geometry we also need the relative positions 
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Figure 4. Graph of the non-negative roots (x = 0, x = x m [ n (Q) and x = £ max (9)) of 
the potential V(x) together with the horizon x = of the "out" Reifiner-Nordstrom 
spacetime. Bounded trajectories are delimited by the x — line and the x = x m i n (Q) 
curve, so that they always cross the horizon of the Reifiner-Nordstrom spacetime, which 
the graph shows to be always smaller than x = x m i n (Q). 



of the horizons with respect to the classical trajectories. This is also briefly discussed 
later on, but it is useful to report here a general result that can be deduced from 
figure |H where for < < O cr it. the turning points, 0, x min , x max and the horizon 
of the exterior metric = G) are plotted as functions of O. It can be seen 

that the classical bounded trajectory, corresponding to the region between the x 
axis and the x m i n dashed curve, crosses for all values of 9 in the considered range 
the exterior horizon at x = 0. The horizon of the internal de Sitter domain (which 
is not plotted and corresponds to the horizontal line x = 1 in rescaled variables) is 
instead never crossed by a bounded trajectory^. 

Asymptotic behaviour : 

quite generally, from (J2Uj) we also see that 



Regularity at the origin : 

the first and second derivatives of the potential are vanishing at x = 0, dV(x)/dx = 
d 2 V(x)/dx 2 = 0, and the third derivative is positive, d 3 V(x)/dx 3 = 6/|9|, so that 
x = is a local maximum for V(x). 

% Nevertheless it equals x max for 8 = 1/2. This is an interesting limiting situation in the case of 
tunnelling across the potential barrier, which will be discussed elsewhere. 



x— »+oo 




(22) 
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Figure 5. The graphical method to study the classical shell trajectories consists 
in plotting the potential V(x), together with the curves for the metric functions of 
the interior and of the exterior. The light gray path is a classically allowed bounded 
trajectory: as discussed in the text it crosses the horizon of the external geometry but 
no changes in the orientation of the normals occur along it. 



Thanks to the above properties, we can now perform the study of the classical dynamics 
for bounded trajectories by restricting the parameter to the range 

3X3/2 



0<|e|<(-j , (23) 

where we have the general situation shown in figure El The figure shows, for = 0.55, 
the potential V(x) together with f- m (x) and f ut(x). The classically allowed path is the 
thicker light gray segment on the x-axis. Let us consider the dynamics in the de Sitter 
spacetime: the shell expands from vanishing radius up to a maximum (grayed path in 
the figure), which remains inside the de Sitter cosmological horizon, since, as can be 
seen, it is not crossed by the trajectory; then the shell shrinks back to zero radius. 
The sign of <7i n does not change along the trajectory, since as we can see always from 
figure there are no points on the trajectory in which the plot of fi n (x) is tangent to 
V(x): moreover, as can be easily verified, it is always positive, so that the trajectory 
crosses the left part of the de Sitter Penrose diagram. This is shown in figure |B1 where 
the interior region is the shaded area, since for a- in = +1 the exterior normal is pointing 
to the right. In the same way we can perform the analysis in the Reifiner- Nordstrom 
domain: we can see (in figure but also from the above discussion about figure EJ) that 
the bounded trajectory during the expansion from a vanishing radius, as well as during 
the following collapse, crosses the horizon of the exterior geometry. As before the sign 
of o~ out does not change along the trajectory (in the zoomed region of figure El we can 
more clearly see that the metric function graph is not tangent to the potential) and 
thus cr out = — 1 always. If we draw the associated Penrose diagram, the exterior region 
is the shaded one in figure [U The results shown above rest on two hypotheses, that we 
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Figure 6. Penrose diagram of the de Sitter interior geometry with the bubble 
trajectory. The interior domain is the shaded region in the diagram. 




Figure 7. Penrose diagram of the Reifiner-Nordstrom exterior geometry with the 
bubble trajectory. The exterior domain is the shaded region in the diagram. 

implicitly took for granted but which deserve a more detailed treatment. 

The first one concerns the stability of the expanding and recollapsing shell against 
single particle decay: if the shell were not stable at the moment of time symmetry, then 
it would become thick and its trajectory would not be approximated by a sharp line (as 
depicted in figures El or |7j) . It is possible to show, following the treatment of jH] (please 
see |Appendix B| for details), that configurations stable against single particle decay of 
charged and/or uncharged particles actually exist. 
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Figure 8. Penrose diagram of the de Sitter interior and the Reifiner-Nordstrom 
exterior geometries joined along the bubble trajectory. It is obtained by joining the 
shaded regions of the two previous figures. The classical dynamics of this compound 
spacetime is described by Israel's junction conditions and is discussed in detail in the 
text. 

The second issue is about possible changes of signs in cr in or <r out , which would 
require a different analysis with respect to the one performed above. We also devote an 
appendix QAppendix Q to show that bounded trajectories are not affected by changes 
of sign in a in or cr out , so that the analysis performed above in a particular case is indeed 
valid in general. 

With these remarks in mind, the complete spacetime manifold can now be 
confidently obtained by joining the interior with the exterior along the shell trajectory, 
i.e. joining the two shaded regions in figure El and figure [7| to get the final result shown 
in figure |H1 An observer inside the shell detects a non- vanishing cosmological constant. 
He lives as an observer inside a cosmological horizon. But as soon as he crosses the shell 
trajectory, he experiences a completely different situations: the cosmological constant 
suddenly vanishes and he can now detect a non-vanishing electric and gravitational 
field, as if outside a body with mass M and charge Q (with M = \Q\, because of 
our simplifying assumptions). We are now interested in studying the properties of this 
gravitational configuration when the system can be considered to be in a semi-classical 
quantum regime. For this we need an evaluation of the classical action along the classical 
trajectory of the shell. 
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4. Numerical evaluation of the classical action 

As already anticipated in (j!2)l we will evaluate the classical action as the integral of the 
classical momentum along a classically allowed trajectory. In our case, remembering 
the naming conventions about the zeros of V(x), this implies that the relevant turning 
points for the bounded trajectory in rescaled coordinates are Ri/H = xi = and 
R-z/H = %2 = x min , so that 



S = A = 2 / P(x)dx 



Sign{(l-x 2 )(x-|e|) 2 } 



The above integral has to be computed when the turning point x min actually exists, i.e. 
in the range for specified by (J2BJ)- This means that the Sign at the exponent is always 
+1, and we can forget about it, so that the above turns into 



* = 2jf P {x)ix ^l xtanh-'j ^^W^ j^ (24) 

We note that when x — 0, i.e. the shell is crossing a (double) zero of the external 
extremal ReiBner-Nordstrom spacetime, the momentum is ill defined, since the argument 
of the inverse hyperbolic tangent is —1. x = G is thus a singularity on the integration 
path, but, being of the logarithmic type, it is integrable (the leading contribution to the 
singularity is determined in Appendix D| ). 



The integral in (|24j) is not exactly computable analytically, but being reassured 
by the considerations above about its existence, the integration can be performed 
numerically. We have performed this kind of analysis with Mathematica® , evaluating 
the integral numerically for 10000 equally spaced test values of O in the interval 
[10~ 4 , (3/4) 3 / 2 ] and for other 10000 test values in the interval [10~ 16 , 10" 4 ] taken as 
a sequence converging to + as l/n A for n — > +oo. The final result is plotted in figure El 



5. WKB Quantum States 

We now assume that the system is in a quantum regime; we will perform its semiclassical 
quantization a la Bohr-Sommerfeld, i.e. considering the action as an integer multiple 
of + h. Remebering that all the computations of the previous section are in terms of the 
adimensional variables defined in (|14|) and (|15|). we have 

S (Q, H) = H 2 S (6) = H 2 S (%) , (25) 



so that we can rewrite the Bohr-Sommerfeld quantization condition as 

S(Q,H) = l\n = n , n = 0,1,2,.... (26) 



+ We work in units where l 2 p = h = c= G= l. 
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5(0) 




Figure 9. Graph of the functional dependence of the action from the parameter as 
results from the numerical integration of l|24|l for 20000 test values of O in the interval 
[10~ 16 , (3/4) 3/2 ] chosen as described in the text. 



5.1. Preliminary Estimate 

We first note that the quantization can be interpreted as giving a relation among the 
parameters of the model, in our case the charge Q and the de Sitter cosmological horizon 
H. Let us take for example a total action of the order of the quantum, /p, i.e. associated 
to a state with quantum number n — 1, for a de Sitter cosmological horizon H = \/2. 
Then S ~ 1/2, i.e. = Q/H pz 0.55, so that Q ~ 0.78. This shows, as a preliminary 
estimate, that a small gravitational system with quantum properties is conceivable. 

We can also see how the action behaves for variations of the parameters Q and H 
by a numerical plot of the level curves of the action. This is shown in figure E3 To 
get a clearer plot in the region close to the limiting line H = (4/3) 3 / 2 Q for small Q, a 
smaller region of the plot is shown, enlarged, in figure ^2 



5.2. Approximating the action 

To get some analytical result we try to fit the 20000 points for which we evaluated the 
action with some simple (polynomial) function. In this way we will be able to get an 
(approximated) relation among Q, H and n. The choice of the approximating function 
is done with two goals in mind: 

(i) to approximate in the best possible way the behaviour of the action S at least in 
some regime; 

(ii) to have a simple enough expression to get an algebraic relation among Q, H and n. 

To fulfill the above requirements we first analyze the leading O dependence of the action 
S, expressed as the integral (j24j) . From |Appendix E we se that P(x) ~ x 2 + O(x) 2 and 
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Figure 10. Graph of the level curves for the action for the level values 1,2,..., 50. 
The levels are obtained using the numerically evaluated action with the function 
ContourPlot of Mathematica®. The thicker dashed black line displays the limit 
H = (4/3) 3 / 2 Q of the condition Q/H > (3/4) 3 / 2 for which bounded trajectories exist. 
The grayed region is shown, blew up, in figure ITT1 
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Figure 11. Blow up of the region of small Q, H parameters. This better shows that 
the quantization condition implies a non-vanishing minimum allowed value for both 
the charge Q and the de Sitter horizon H , i.e. a non- vanishing minimum allowed value 
for the charge Q together with a maximum allowed value for the cosmological constant 
l/H. 
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Figure 12. Comparison between approximated and numerical evaluated actions. 
The expression 5 rog .(0) is obtained, as described in the text, by a fit to a suitable 
polynomial. This figure aims to show that the main behaviour can be qualitatively 
approximated, but must not be taken to imply that the approximation is everywhere 
very good. Indeed when the action is very small, a very small approximation error can 
still give a relevant relative uncertainty. 



the upper turning point x min (0) ~ O + 0(9), so that S{Q) ~ 3 + O(0) 3 . We thus 
choose an approximating function starting with a third power of and having the next 
two powers of 0: 

Ax 5 + Bx A + Cx 3 , (27) 

The coefficients A, B, C in ()27j) have been determined from the 20000 sample points 
using the function Regress of Mathematica®, with the parameter IncludeConstant 
— > False, since there is no constant term in the model: they result to be, together with 
the corresponding standard errors, 

A = +10.11 ±0.04 

B = -7.39 ± 0.04 (28) 
C = +3.64 ±0.01, 

with an adjusted regression coefficient of 0.99981. 

The comparison between approximated and numerical evaluated actions is plotted in 
figure Level curves of the approximated action are plotted in figure 

From the approximated expression ()27j) , which by ({251) anc ^ mus t equal n when 
multiplied by H 2 , we can get the following equation of the third degree in H, 

AQ 5 + BQ 4 H + CQ 3 H 2 - nH 3 = : (29) 

this can be solved exactly to get the approximated relation for if as a function of Q 
and n as it comes from the simplified expression (j2T|) for the action. The quite complex 
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Figure 13. Graph of the level curves for the approximated action levels 1,2,..., 50. 
The levels are obtained using the approximating polynomial action with the function 
ContourPlot of Mathematica®. The thicker dashed black line displays the limit 
H = (4/3) 3 / 2 Q of the condition Q/H > (3/4) 3 / 2 for which bounded trajectories exist. 



algebraic expression is 
H(Q;n) = — + 



3nVH 2 (Q; n) + y/4Q*H*{Q; n) + H*{Q; 



n 



Q 5/3 VH 2 (Q;n) + y/4QW*(Q;n) + H*(Q; 



n) 



2 1 /33 n 

where 

Hx(Q;n) = -C 2 Q 2 - 3B& 

H 2 (Q; n) = -2C 3 Q* - 9BCnQ 2 - 27 A 2 n 2 

and the approximated relations above are only valid if 

101 



H > 



(3/4) 3 / 2 : 



(30) 



a condition plainly coming from (|23J) . The approximated levels of figure EI] are nothing 
but the graphs of the n = const, relations coming from ()3()j) . 

6. Discussion 



We have presented a model in which it is analyzed a general relativistic system composed 
of two spacetime domains, a de Sitter interior with cosmological horizon H and an 
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extremal Reifiner-Nordstrom exterior with M = \Q\, joined across a thin shell of mass- 
energy m — \Q\. A semiclassical quantization of classically bounded trajectories can 
be performed using a scheme a la Bohr-Sommerfeld. In this way it is found that the 
quantum states can be characterized by a quantum number n, which is responsible 
for restricting the allowed values of Q and H. In particular the quantum dynamics 
of the system constrains the values of H as functions H = H(Q\ri), which we have 
approximated after a numerical analysis of the problem. 

It is interesting to note that this is a sensible result in a semiclassical approximation 
to a quantum gravitational situations. Indeed a full quantum theory of gravity would 
treat the three-metric as a dynamical infinite-dimensional variable to be determined, 
say, by the Wheeler-deWitt equation in superspace. From our point of view, under 
the assumption we made, we are only in a minisuperspace approximation, since the 
functional form of the metric functions is fixed from the very beginning, leaving H and 
Q as the only free parameters: quantum gravity will thus impose condition on these 
quantities, i.e. determine the only residual degrees of freedom in the three geometry. 
We can also see from the enlarged plot of figure El that the quantization condition 
selects a minimum value for the allowed charge of the semiclassical quantized shell. 
At the same time, the cosmological constant of the interior region, cannot exceed a 
maximum value, which is also a sensible consequence of quantization. For a given fixed 
value of the ratio Q/H = k, which means that we are "moving" across the graphs 
of figures El and El along a line through the origin, we have discrete allowed values 
for Q and H: these point are located on a parabola, since from (|25jl we exactly have 



We also note that by the final result, systems characterized by different scales in 
the parameters can be described: among these, as expected, we find small scale systems 
(with a charge which is a small multiple of the elementary electron charge). Moreover 
an external asymptotic observer measures a total mass energy for the shell given by 



outside domain an object with rest mass m = \Q\ and charge Q, which is the exterior 
manifestation of a bounded interior containing a part of spacetime characterized by a 
non- vanishing vacuum energy. Due to the form of the potential this bound semiclassical 
state is metastable, and will decay into an infinitely expanding shell after a finite time 
(a similar situation occurs in |%T|). which in principle could be calculated (as proper 
time), studying the process of tunnelling across the classical effective potential barrier*. 

* This will be the topic of a forthcoming paper. 



S(iiH, H) = H 2 



S(k). 




which in our case is nothing but E — m, since M 2 



Q 



2 : thus we effectively see in the 
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Appendix A. Zeros of the potential and critical value of 

The zeros of the potential V(x) can be obtained in closed form, since they are the zeroes 
of the numerator, i.e. the solutions of the equation 

V(x) = => x 2 (x 4 -Aax + Aa 2 ) =0 (A.l) 

with a > and — x 3 + 2x — 2a ^ 0. We are interested in the non-negative solutions, 
which, apart from the x = one, can be determined exactly as solutions of the fourth 
order equation 

x 4 - Aax + 4a 2 = 0. (A.2) 

By setting 

^=| 9 a 2 + (81a 4 - 192a 6 ) 1/2 } 1/S 
we have that for < a < (3/4) 3/2 



1 I /4-3 1 /3 a 2 



4 ■ 3V 3 fl 2 6 • 2 1 /2 r 

B+ , =}. (A.3) 



3 J :S 2 J 2 I V B \\ a ^/ 4.3V3 a2 | B 

These expressions are not very enlightening: the one for x min has been used to exactly 
evaluate the upper integration limit in the numerical evaluation of the integral that 
gives the classical action. 

To see when the quartic part of the potential has two positive roots we can use the 
expression above, but also a smarter procedure, as follows. Clearly the limiting case is 
the one in which the potential is tangent to the positive x axis, i.e. the two solutions 
coincide. In this case the quartic part must be of the form 

(x - a) 2 (x 2 + fix + 7) 

= x A + (J3 - 2a)x 3 + (7 - 2af3 + a 2 )x 2 + (a(3 2 - 2a^)x + a 2 7, 

which by comparison with x A — 4ax — 4a 2 gives the set of equations 
P - 2a = 
7 - 2a(3 + a 2 = 
a 2 (3 — 2«7 = —4a 
a 2 7 = 4a 2 

Then a = (3/2 and 7 = 3/? 2 /4 from the first two equations. This gives a = a 1//3 from 
the third and a = (4a 2 /3) 1//4 from the fourth. These last two relations are compatible 
for non-vanishing a, if and only if a = (3/4) 3 / 2 , which is thus the critical value of a. 



Appendix B. Stability of the shell against single particle decay 

In this section we discuss the stability of the trajectory of the infinitesimally thin shell 
under single particle decay, following the treatment that can be found in jH]. The 
proof that the shell is stable against single particle decay of uncharged particles is not 
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reproduced here, because it can be easily derived from the reference cited above, to 
which the reader is referred. We will instead shortly discuss the case in which charged 
particles are involved. 

In more detail, the relevant question is if the motion of a charged particle, which at 
the instant of maximum expansion starts out where the shell is located, will subsequently 
be governed by a confining, i.e. "U-shaped", effective potential, or not. Performing the 
analysis at the instant of maximum expansion, where the potential is static, simplifies 
the computation: subsequent changes of the potential will have, anyway, only adiabatic 
effects on the locally trapped particle and this is not relevant for the point under 
discussion. To get the desired result we will proceed in two steps: 

(i) identify the effective potential governing the motion of a particle in the exterior 
ReiBner-Nordstrom geometry; 

(ii) evaluate if it is "U-" or "fl-shaped" at the point of maximum expansion. 
We will work in the adimensional units used throughout the rest of the paper. 

Appendix B.l. Effective potential for a charged particle in the Reifiner- Nordstrom 



The effective potential for the motion of a particle of charge q and mass fj, in the 
ReiBner-Nordstrom spacetime can be obtained in many different ways. Probably the 
quickest one is to start from the more general result that can be found in the equation 
after equation (3) in Box 33.5 of [70], i.e. the effective potential for the orbits of test 
particles in the equatorial plane of a Kerr-Newman black hole. Specializing this result 
to a black hole with zero angular momentum we obtain the effective potential in the 
ReiBner-Nordstrom case. Perhaps more instructive is to perform again the analysis until 
equation (6) of [B] adding the electrostatic contribution to the particle momentum or 
solving the Hamilton- Jacobi equation in the ReiBner-Nordstrom metric. Anyway, the 
final result for the extremal case we are interested in is 



Appendix B.2. Evaluation of the Effective Potential at the point of maximum expansion 

With the above results at hand, we now evaluate the second derivative of the effective 
potential, i.e. 



spacetime 





where e is the charge/mass ratio of the test particle and A = L/H is its angular 
momentum per unit mass (L = L/fi) in the H scale defined in (1 14 J) and ([TKjl . 



V(x) 




(B.2) 



3/2 
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we are interested in its sign at the value x m i n , the maximum radius of the shell given in 
flA.3|) : a positive sign will indicate that the potential is "U-shaped" and thus the shell 
stable. To get the result, we can restrict the study to the numerator iV(£ m i n ) of (|B.2|) . 

/ _ \ 3/2 

N(x min ) = 3x 2 min \ 2 + 2A 4 + 2e0 (x^ + A 2 J - (2x 4 min + 9x 2 min \ 2 + 6A 4 )|0|, 

since the denominator does not contribute to the sign. We also remember that the 
quantity x m i n depends on 9, since (|A.1J) comes from (|2T)|) with a = |0|. 

Even with the above simplification, the detailed study of the sign of the quantity 
under consideration is complicated, mainly because of the non-trivial dependence; 
a graphical analysis is also of little help, since N(x m i n ) depends on the three variables 
0, e and A, namely the adimensional charge Q/H, the charge/mass ratio q/fi of the 
test particle and the adimensional angular momentum per unit mass L/(Hfi) of the 
test particle. Thus we will not search for the most general result, i.e. we will not give 
necessary and sufficient conditions for the stability of the shell; we will show, instead, 
that in some physically reasonable situations the shell itself is indeed stable against 
single charged particle decay. 

In particular we see that in the units we are using, the adimensional parameter e 
is a large number, i.e. the charge/mass ratio for an elementary particle is very large. 
Thus we consider the behaviour of the numerator for large e: 

lim N(x min ) = +oo; 

e— >oo 

the limit has always the "+" sign, since we consider emission of charges of the same 
sign of those composing the shell; thus in this limit the second derivative of the effective 
potential is positive, which shows stability under charged elementary particle decay. 
As a second case, we see what happens for radial emission of particles: 

lim N(x min ) = 2a4 in |0|(|e| - x min ); 

we again used the fact that ejected particles have the same charge as the shell, so that 
e0 = |e| • |0|. In this case also we see that, for elementary particle emission, we certainly 
can realize the situation e > x min , so the sign is again positive. 

Even restricting the study to the two cases above, we can thus conclude that, shell 
configurations which are stable against single particle decay can be realized. 

Appendix C. General analysis for a in , a ont on a bounded trajectory 

The "graphical" analysis of the classical dynamics performed in section E] is based on 
the plot of figure El In the discussion a relevant point is the sign of <7 in and cr out , which is 
crucial in determining the direction of the normal to the shell pointing in the direction 
of increasing radius as well as the side of the Penrose diagram crossed by the trajectory. 
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We will show here that the situation analyzed for the particular value = 0.55 is, in 
fact, general. The sign of a in is given (fTSj). so that we see that is positive for 

x < x^{Q) = {2Q) l l\ 

More complicate is the analysis of the sign of a ont , because from (fT^j) we see that it is 
determined as the sign of a polynomial of order four. It has at most two real roots for 
< ©crit./v^ = (3/4) 3 / 2 /a/2; this can be seen also from the analysis of |Appendix A 
Indeed a out > if 

x 4 - 2ax + 2a 2 < with a > 0. 

The real roots of the left hand side of the above inequality can be deduced from those 
of the left hand side of ()A.2jl since 

x — 2ax + 2a — ► - (x — Aax + 4a when x — > — = and a — > — =. 

4 v ' ^2 V2 

Then, if we call x^~l t and x^2 t the two real values for which <j out changes sign, we see, 
using the transformation above, that we must have a < (3/4) 3//2 /v^2 and that 



with 



8 ' 3W T + ; 125 > (C.l) 



\ 7 ^p + T 



T = y 186 2 + 2V8 16 4 - 38 46 6 

Since the above expressions are not very enlightening (and the same is true for those of 
(|A.3|l ). the easiest way to compare them with the turning points is again the graphical 
one, i.e. the plot of x CTin (9), x^ u \(6), x min (6), x max (9), which we can see in figure IUT1 
For small x (on the vertical axis) <r in is positive and a out is negative. Since all the zeroes, 
when they exist, are bigger than the the turning point x min , which is the upper limit of 
the bounded trajectory, then there is no change of sign of cr's along it. 

Appendix D. Character of the singularity on the integration path 

In this section we determine the leading contribution to the logarithmic (and thus 
integrable) singularity on the integration path that appears in the evaluation of the 
integral in (|24|). The logarithmic character stems from the definition of the inverse 
hyperbolic tangent in terms of the logarithm and from the fact that its argument is a 
rational function with the following properties}]: 



r ,, .. ny\x;Q) x^ 4 -4|0|x + 46 2 
hm J-{x: = hm — — — — ^— = hm ; — — = 1, 

x->\e\ v ' ; x-|e| V(x;Q) s-[e| x 3 -2x + 2|0| 

i.e. 

^(|e|;6) = l (D.l) 
jj We define T , 72. and V according to the first two ='s of the equation below. 
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Figure CI. Graph of the curves that helps in determining the sign of a's for different 
values of the parameter 0. Since the classically allowed trajectories are delimited by 
the x = axis and by the dark gray continuous curve x = a; m i n .(0), we see that no 
changes of signs of sigma's occur along them. 



To extract the behaviour of the above function around the point x — |9| we develop it 



in power series. Let us set 

K{x;Q) =x 6 -4|9|:r 3 + 4eV with ft(|e|;9) = e 6 (D.2) 

V(x;Q) =x 3 -2x + 2\&\ with P(|6|; 0) = |0| 3 . (D.3) 
Then it follows 

n'(x;B) = Qx 5 -I2\&\x 2 + 8& 2 x so that K'(\Q\; 6) = 6|0| 5 - 4|0| 3 (D.4) 

V'(x;&) = 3x 2 -2 so that £>'(|9|; 6) = 30 2 - 2. (D.5) 

and 

n"(x;Q) = 3Ox 4 -24|0|a; + 80 2 so that K"{\Q\; 0) = 3O0 4 - 160 2 (D.6) 

V"(x; 0) = Qx so that £>"(|0|; 0) = 6|0|. (D.7) 



We then compute JF'(|0|; 0) and JF"(|0|; 0), i.e. the first and second derivatives of the 
argument of the inverse hyperbolic tangent. Since 
_ K'V - 2TZV 

we obtain 

jr'(|0| ; 0) = o. (D.8) 
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Moreover 



r 



■a 



2 (K"V - K'V - 2TZV") TLV - (K'V - 2KV) (K'V + AKV) 



47^3/2^3 



and we get 



^ ,, (i i5 ) = -4( 1 - 02 )- 



(D.9) 



The expansion of T around x = |6| up to second order can then be written, using (|D.1|) . 
(|DHjl and dEU), as 



and inserting this result inside the expression of the hyperbolic tangent in terms of 
logarithms, we can easily see, as expected, that the singularity is integrable. 

Appendix E. 6 dependence of the Action 

We consider the action integral f|24jl . Expanding the integrand we get 



We thus see that the leading term for small |0| is ~ |© 3 |- 
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